Transverse rectification of disorder-induced fluctuations in a driven system 
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We study numerically the overdamped motion of particles driven in a two dimensional ratchet 
potential. In the proposed design, of the so-called geometrical-ratchet type, the mean velocity of a 
single particle in response to a constant force has a transverse component that can be induced by 
the presence of thermal or other unbiased fluctuations. We find that additional quenched disorder 
can strongly enhance the transverse drift at low temperatures, in spite of reducing the transverse 
mobility. We show that, under general conditions, the rectified transverse velocity of a driven particle 
fluid is equivalent to the response of a one dimensional flashing ratchet working at a drive-dependent 
effective temperature, defined through generalized Einstein relations. 

PACS numbers: 74.25.Qt,05.40.-a,05.45.-a,05.60.Cd 



The idea of generating a directed dissipative transport 
in a system kept out of thermal equilibrium only by un- 
biased perturbations has motivated an outburst of ex- 
perimental and theoretical works in the last years. The 
ratchet effect is indeed of interest, both for applications 
and modelling, in very diverse systems, ranging from 
biological motors, 1 colloidal matter^ granular matter^ 
vortex matter in superconductors^ Josephson junction 
arrays^ atoms in optical traps,.^ electrons in semicon- 
ductor heterostructuresji to gambling gamesi& One of 
the simplest models is the "flashing ratchet" , where a 
directed motion of a Brownian particle (i.e. breaking of 
the detailed balance condition) is obtained by coupling it 
to a pulsating asymmetric-periodic potential. The identi- 
fication of the essential physical ingredients for the effect 
shows that a large variety of ratchets and rectification 
mechanisms can be realized^ 

Recently, there has been a growing interest in the so- 
called geometrical ratchets since they can be used as con- 
tinuous "molecular sieves" to separate particles experi- 
mentally (such as macromolecules or mesoscopic objects), 
according to its physical properties. These devices are 
typically two-dimensional systems containing a periodic 
array of asymmetric obstacles. By driving the particles 
through the array an average lateral drift appears, as 
transverse diffusive motion is rectified by the collisions 
with the asymmetric obstacles. Different types of geo- 
metrical ratchets have been analyzed in the literature, 
both experimentally 10 and theoreticallyiii The effect of 
additional quenched disorder in these two-dimensional 
systems has not been discussed yet, though interesting 
anomalous transport properties of one-dimensional dis- 
ordered ratchet systems were reported^ Such a study is 
not only relevant for applications where disorder can not 
be avoided, but it is also an interesting and challenging 
issue. The driven motion of particles in a disordered sub- 
strate yields a non-trivial hydrodynamics. The current- 
driven motion of vortices in type II superconductors is 
a prominent example, where disorder, apart from reduc- 
ing dissipation, is responsible for marked non-equilibrium 
transport and magnetic properties^ On the other hand, 
already the simpler case of driven non-interacting Brow- 
nian particles in two dimensions displays complex phe- 




FIG. 1: (a) Ratchet potential with disorder, (b) Schematics of 
the transverse rectification mechanism (top view). Particles 
move with an average velocity V y in the direction of the ap- 
plied force F y . In the white regions the interaction with the 
(attractive or repulsive) centers and with the thermal bath 
induces diffusion in the non-driven direction In the shaded 
regions a periodic-asymmetric potential tends to localize par- 
ticles at its minima. An average transverse shift (see circles) 
is produced at a rate V. 



nomena. While diffusion is anomalous at equilibrium, 14 
under a finite drive diffusion becomes normal in the co- 
moving frame, with anisotropic and velocity-dependent 
diffusion constants and mobilities^ Moreover, a disor- 
dered substrate can provide alone a local or global ratchet 
effect, such as the generation of large-scale vorticity in the 
probability current by driving particles with an uniform 
alternate drive^ and the net directed motion produced 
by driving the particles with crossed ac-drives^i 

In this paper we investigate the effect of quenched dis- 
order in a simple geometrical ratchet design, under an 
uniform and constant driving force. We find that disorder 
can strongly enhance the transverse drift both for non- 
interacting and interacting particles, thus improving the 
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performance of the device for applications. We show that 
the transverse velocity of a driven fluid is equivalent to 
the response of a one-dimensional flashing ratchet work- 
ing at a drive-dependent effective temperature, defined 
through generalized fluctuation-dissipation relations. 

Let us consider the overdamped motion of particles in 
a two dimensional potential like the one depicted in Fig^ 
The equation of motion of a particle in position R; is: 



dR, 



= -V,; 



+ F + C<(t), (1) 



where Rij = |R ? ; 



Rj| is the distance between parti- 
cles Rip is the distance between the particle i and 
a site at R p , r\ is the friction, and F = F y y is the 
driving force. The effect of a thermal bath at tem- 
perature T is given by the stochastic force Ci(t)i sa ^ _ 
isfying (Cf(i)) = and (tf(f)tf'(O) = 2r,k B T6(t - 
t')SijSau'j where (...) denotes average over the ensem- 
ble of Ci- For concreteness we consider a logarithmic 
repulsive particle-particle interaction V(r) = —A v ln(r) 
which corresponds for instance to the vortex-vortex 
interaction in 2D thin film superconductors^. Par- 
ticles interact with the quenched potential U (R) = 
C/r(R) + U P (R). U R is a ratchet potential with the form, 
U R (R) = £F R (Y)G R {X), where X = R.x, Y = R.y, 
G R (X) = sm(2nX/a) + 0.25 sin(47rX/a) and F R (y) = 
Uq cos(27Tj//6)G[cos(27rj//6)], with the Heaviside func- 
tion. This ratchet potential is similar to a periodic ar- 
ray of obstacles, asymmetrical around the x — axis but 
symmetrical around the y — axis, as the ones consid- 
ered in Ref.iii Disorder is short-range correlated, and it 
is modeled as a random distribution of centers such that 

U r (Ri) = J2 P V (JWr,)a > where R w = |R*"Rpl is the 
distance between the particle i and a center at R p . Cen- 
ters can be either attractive A p < (wells) or repulsive 
A p > (humps) or a combination of both. We solve EqQ] 
numerically by using the method of Ref^£. Length is nor- 
malized by r p , energy by 2 A p , and time by r = rjr p / \ 2A p | . 
We consider N — 60 particles and N p pinning centers in 
a rectangular box of size L x x L y and periodic boundary 
conditions, with L y = 100, L x — 2Q^/2>L y . We average 
calculated properties over 500 disorder realizations. 

We start by discussing the simplest case of non- 
interacting particles, A v = 0, without disorder, N p = 0, 
and with a ratchet potential of amplitude Uq = 1. The 
dashed-lines of Fig^a) show the transverse drift rate 
v = (jr H 1 ) at T = °- 05 as a function 01 the lon- 
gitudinal velocity V y = (jfJ2i^)- We see that the 
transverse velocity V increases from zero, has a max- 
imun V ~ 0.0075 &tV y ^ 0.5, and decays to zero at large 
longitudinal velocity. Since V = at T = 0, the av- 
erage directed transverse motion observed is induced by 
the thermal noise. This rectification effect is easy to un- 
derstand, and Fig^b) illustrates the mechanism, where 
the transverse diffusion constant is D = 2T if there is no 
disorder nor inter-particle interactions. For our discus- 
sion it is useful to make explicit the connection between 



the type of response shown in Fig|2a) and the one of 
a flashing ratchet. If F y is large, the mean velocity in 



the driven direction V y = (^r) 



iS Vy ~ Fy ~ OiF- 1 ) 



and longitudinal fluctuations are much smaller (by a fac- 
tor 0{F y 1 )) than transverse fluctuations^ At T = 0, 
the equation of motion for the coordinate X of a particle 
located at R = Xx + Yx can be thus written as 



dX 



-F R (V y t)G' R (X). 



(2) 



Since F R (V y t) = Uq cos(27rV y t/b)e[cos{2irV y t/b)], X 
feels the ratchet potential G R (X) switching on and off pe- 
riodically with time periods T on = T off = b/2V y . At small 
drives the mechanism is the same although r on becomes 
increasingly larger than r off since the wells of U R (X, Y) 
delay the motion in the y direction, (see FigQ. The 
mapping to a flashing ratchet explains the observed di- 
rected transverse motion with V > when T > and 
can be thus used as an effective model to explain all the 
features of the response shown in FigEfa). 

The effect discussed so far is similar to the one de- 
scribed in Ref.»ii Let us now add disorder, by putting 
N p = 2000 randomly located pinning sites. The resulting 
response is shown in Fig|2Ja) (symbols). As we can see, 
disorder strongly enhances the rectification at intermedi- 
ate and large longitudinal velocity and also broaden the 
range of V y where response is appreciable, with respect to 
the clean case. In addition, we find that at intermediate 
and large Vy, V is finite even in the T = limit, since dis- 
order induces transverse diffusion when V y > 0, even in 
the absence of thermal fluctuations. Finally, let us now 
turn-on the repulsive interaction between particles. In 
FigEtb) we show V as a function of V y , for different val- 
ues of the repulsion strength A v . In the inset of Fig|2f b) 
we see that the maximun response V max is almost con- 
stant with A v up to values A v ~ 0.2 where a slow decay 
starts, but it is larger than the response of the clean sys- 
tem up to A v = 2. The decay of the response at large A v 
is explained by the decrease of transverse wandering due 
to increasingly correlated collective motion^ In Figl^b) 
we show that the response for purely attractive pinning 



centers A v 



T is smaller than for repulsive centers 



A p > 1 for small values of V y , but indistinguishable for 
larger values of Vy. This is due to the fact that, at the 
density of centers considered, attractive centers are more 
effective to pin particles than humps, since the latter can 
provide two-dimensional pinning only by forming rare ge- 
ometric traps. However at a density N p /L x L y ~ 1/r? all 
these differences disappear completely. 

In order to understand the rectification characteristics 
described above it is instructive to study, separately, the 
motion of particles in the purely disordered case, with- 
out the ratchet potential (i.e. Uo — 0). For simplicity 
we consider only the case of non- interacting particles, 
but we expect similar results for interacting particles 
in the dynamical regimes where transverse diffusion is 
non-zero^*2S We analyze in detail the non-equilibrium 
transverse fluctuations as a function of F y , since they 



3 



> 0.010 




> 0.010 

0.005 
0.000 



0.15 



0.10 



0.05 



(a) 




0.3 



0.2 
o 



FIG. 2: Transverse velocity V vs longitudinal velocity V y . (a) 
Symbols correspond to the disordered system at T = 0.05, the 
dashed-line to the clean system at T = 0.05 and the solid-line 
to the clean system at an effective temperature T eii (Vy,T). 
(b) V vs V y for different interaction strengths A v : A v = 0.05 
(+), A v = 0.2 (*), A v = 0.5 (o), A v = 1 (A), A v = 2 (□), 
A v — 5 (x). (•) symbols correspond to purely attractive pin- 
ning centers, with A v — —0.5, at T = 0.05. Inset: maximun 
rectification as a function of A v . 



affect directly the rectification in the presence of the 
ratchet potential. Following Refii& we define the observ- 
ables 0(t) = £ Y£i SiXi{t) and 6{t) = «<*<(*), 
where Si = — 1,1 are random numbers with Sj = 
and SiSj = 5ij. The quadratic mean displacement can 

be written as A(t,t ) = ^ E^'id^iW ~ ^fa)l 2 ) = 

C(t, t)+C(t ,t )-2C{t, t ), with C(t, t ) = (O(i)O(to)). 
The integrated response function x for the observable O 
is obtained by applying a perturbative force fi = eSji; at 
time to and keeping it constant for all subsequent times 

l 



on each particle, x(Mo) = lim e _> 7 (O(t)) e -(O(t)) e=0 
In the steady state A(i, to) — A(t — t ) and x(t, f ) = 
x(t — to) and in particular at equilibrium the fluctuation- 
dissipation theorem (FDT) imposes x(t) = A(t)/2T. 
When F y > the system is out of equilibrium and the 
FDT does not hold. We will show however, that gen- 
eralized fluctuation-dissipation relations can still be de- 
fined for our system. In the long time limit we find 
A(t) ~ Dt and x(t) ~ ^t thus allowing us to define the 
transverse diffusion constant D and the transverse mo- 
bility [i. These two quantities depend on the longitudinal 
driving force as shown in Fig[3{b). D is non-monotonic, 
has a peak at F y ~ 1.5, and decays approximately as a 
power law towards the equilibrium value without disor- 
der 2T for large forces. This behavior can be understood 
by considering the effective transverse random walk in- 
duced only by the collisions with the pinning centers at 



FIG. 3: Motion in a purely disordered potential as function 
of F y . (a) Transverse (x direction) diffusion constant D and 
mobility /j,. Dashed and solid lines indicate asymptotic forms 
of D. (b) Effective temperature T ett ~ D/2/i. The dashed- 
dotted line indicates the thermal bath temperature, T = 0.05. 
The inset shows that T m u satisfies a generalized fluctuation- 
dissipation relation between the integrated response \{t) and 
the quadratic mean displacement A(i). The dashed line cor- 
responds to the equilibrium fluctuation-dissipation relation, 
which is valid only at short times t < (,/V y . Upper symbols 
correspond to F y — 0.0, and lower symbols to F y = 4.0. 



T = 0, and by simple heuristic arguments is possible to 
find the asymptotic forms D ~ n v r^V y at small V y and 
D ~ n p rpAp/Vy at large V y r^ indicated in FigEJb). At 
large F y the transverse mobility fi approaches the equilib- 
rium value without disorder, /j, = 1 (independent of T). 
At small F v , /i decreases due to trapping and its value 
at the limit F y — >■ is controlled by T. In the inset of 
FigEH c ) we show the parametric plot of x(t) vs A(t)/2T 
for Fy = (equilibrium) and F y = 4.0 (out of equilib- 
rium). We see that the equilibrium FDT holds for F y = 
as expected. For F y = 4.0 (and in general for F y > 0) we 
observe instead that the FDT holds only at very short 
time scales, t < r p /V y , but it is violated at long times. 
The type of violation observed can be quantified using 
the notion of time-scale dependent "effective tempera- 
tures" introduced by Cugliandolo, Kurchan and Peliti. 21 
Following Ref:— we define a velocity-dependent trans- 
verse effective temperature T eff from the slope shown in 
the inset of Fig|3c) . At long times this implies the gen- 
eralized Einstein relation T eff ~ D/2/i in the non-driven 
direction. In Fig|3|[c) we see that T eff follows closely D 
except at low forces where \i decreases towards the value 
fi ~ 0.5 at very low forces. For interacting particles simi- 
lar results for T eff were obtained at low and intermediate 
forces, in the plastic and smectic regimes of motion^ At 
large forces and small temperatures however, the forma- 
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tion of pinned rough channels for particle motion^ leads 
to a transverse freezing of the moving vortex fluids at 
T = 0. At this dynamical transition D and \i are strongly 
reduced, and vanish at T — strictlyii^i^ 

The fluctuating motion in a purely disordered sub- 
strate and the connection between the response of a geo- 
metrical and a flashing ratchet discussed above motivate 
a simple model that can be used to describe the rectifica- 
tion response of the disordered geometrical ratchet. We 
propose the following equation for the transverse motion 
of non- interacting particles, 

1 riX 

^^n-FntyMm + ^it). (3) 

where we have replaced Yi ~ V v t as before, the pin- 
ning force by an effective thermal noise Q ff (t) and the 
bare mobility by fJ>(V y ). We use that (Qfi(t)} = 
and (Q lt (t)& f (0) = (2T eli (V y )/ l i(V y ))6(t-t')S ij where 
T e a (T, V y ) and ^{T, V y ) are the ones shown in Fig|3fa). 
Eq. therefore models transverse motion in a coarse- 
grained way (in time and space), satisfying the gener- 
alized fluctuation-dissipation relation shown in the inset 
of Fig|3Ib) (in the absence of the ratchet potential). By 
construction, the assumptions of the model are that (i) 
transverse forces are small compared with the longitudi- 
nal drive F y and (ii), the particle motion is incoherent 
at the length scales of the ratchet potential. In Figl^a) 
we see that the transverse drift generated by this model 
is close to the one of the full model, for the parameters 



analyzed in this paper. The rectification characteristics 
of the two-dimensional geometrical ratchet are therefore 
well described by the one-dimensional flashing-ratchet 
described by Eq[3] working at the effective temperature 
7eff(T, V y ) and friction /i _1 (T, V y ), determined by the 
disorder and the longitudinal velocity. Using this model 
the enhancement of the rectification observed in Fig[2la) 
can be simply attributed to the fact that T eff > T, but 
with T eff still smaller than the optimal temperature for 
rectification of the effective pulsating ratchet Eq|2|is also 
expected to work for interacting particles by using the 
respective T eff (T, V v ) and /i(T, V v ) except at low T and 
large F y where condition (ii) can be violated since parti- 
cles become correlated over long times and distances^ 

In conclusion, we have studied numerically the effect 
of quenched disorder in a geometrical ratchet. We find 
that disorder enhances the transverse rectified velocity 
of a driven fluid. If particle motion is incoherent at the 
scale of the ratchet potential the response can be simply 
described by a one-dimensional flashing ratchet working 
at a disorder-induced, drive-dependent effective mobil- 
ity and temperature, satisfying generalized Einstein rela- 
tions. This effect can be used experimentally to enhance 
and control the performance of geometrical ratchets. 
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